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ABSTRACT 

We report the discovery of the correlated optical/X-ray low-frequency quasi-periodic 
oscillations (QPOs) in black hole binary SWIFT J1753.5-0127. The phase lag between 
two light-curves at the QPO frequency is close to zero. This result puts strong con¬ 
straints on the nature of the optical emission in this object and on the origin of the 
QPOs in general. We demonstrate that the QPO signal and the broadband variability 
can be explained in terms of the hot accretion flow radiating in both optical and X-ray 
bands. In this model, the QPO appears due to the Lense-Thirring precession of entire 
flow, while the broadband variability in the optical is produced by two components: 
the hot flow and the irradiated disc. Using the phase-lag spectra, we put a lower limit 
on the orbital inclination i > 50°, which can be used to constrain the mass of the 
compact object. 

Key words: accretion, accretion discs - black hole physics - stars: individual: SWIFT 
J1753.5-0127 - X-rays: binaries. 


1 INTRODUCTION 

Accreting black holes (BHs) in our Galaxy have been at a 
close look by the X-ray telescopes for the last half a century. 
Most of them are transient sources, typically being active 
for a few months and spending years to decades in qui¬ 
escence (Remillard & McClintock 2006). X-ray studies re¬ 
vealed that the objects go through several spectral states 
during the outburst, classified according to the energy range 
where most of the energy is liberated (Zdziarski & Gierlihski 
2004) . BH X-ray binaries demonstrate a tight connection be¬ 
tween the spectral states and the variability properties (see 
e.g. Gilfanov et al. 1999; Axelsson et al. 2005). The objects 
are highly variable on the sub-second time-scales in the hard 
state and this variability is largely suppressed when it enters 
the soft state (see review hr Done et al. 2007). 

The accretion geometry close to the compact objects 
and physical mechanisms shaping their broadband spectra 
are among the major questions still being debated. Analy¬ 
sis of X-ray spectral and timing properties are the power¬ 
ful tools to probe the physical processes operating in the 
vicinity of compact object. The shape and the magnitude 
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of the power spectra are connected to the spectral states 
and, ultimately, to the accretion geometry. In the hard and 
hard-intermediate states, the variability spectra of BHs are 
characterised by a high amplitude broadband noise (rms 
amplitude of about 30 per cent), in addition to which nar¬ 
row features known as quasi-periodic oscillations (QPOs) 
are often detected (Belloni & Stella 2014). Hereafter, the 
QPOs refer to the low-frequency QPOs observed in the fre¬ 
quency range ~0.1-10 Hz, with the central frequency corre¬ 
lated to the low-frequency break of the broadband noise and 
the X-ray flux (Wijnands & van der Klis 1999; Casella et al. 
2005). The whole set of the QPO models was developed, yet 
the connection of these timing features to spectral forma¬ 
tion has not been addressed in most of them. The most 
established in this respect is the model where QPOs ap¬ 
pear due to the Lense-Thirring precession of a hot accretion 
flow (Stella & Vietri 1998; Fragile et al. 2007; Ingram et al. 
2009). 

Lately, intriguing properties of the fast optical variabil¬ 
ity and its correlation with the X-rays in Galactic BHs with 
low-mass companions also captured attention of the commu¬ 
nity. This variability is likely produced in the processes of 
accretion/ejection and thus can also serve as a probe of the 
accretion physics. The most straightforward method to find 
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the connection between optical and X-rays is to calculate the 
cross-correlation function (CCF). Three transients, GX 339- 
4, XTE J1118-f480 and SWIFT J1753.5-0127 were found 
to demonstrate an anti-correlation between light-curves in 
these energy bands (Motch et al. 1983; Kanbach et al. 2001; 
Hynes et al. 2009; Durant et al. 2008; Gandhi et al. 2008, 
the so-called precognition dip in GCF). These observations 
immediately ruled out the reprocessing nature of the opti¬ 
cal photons, favouring their non-thermal origin, in the hot 
accretion flow or, perhaps, the jet (see e.g. recent reviews 
by Poutanen & Veledina 2014; Uttley & Casella 2014). In¬ 
terestingly, many observations also revealed the presence of 
the QPOs in the optical, ultraviolet and infrared power spec¬ 
tra at frequencies similar to those in the X-rays (Hynes et al. 
2003). So far, the relation of the optical broadband variabil¬ 
ity as well as the QPOs to the noise and low-frequency QPOs 
in the X-rays was not studied. 

In this paper we study the variability patterns of the 
BH binary SWIFT J1753.5-0127. We use the data obtained 
in 2007, when the source was in the hard state. The simul¬ 
taneous observations report the presence of the QPO in the 
optical power spectra, but no statistically significant QPO 
in the X-rays was detected (Durant et al. 2009). Using the 
CCF and the phase-lag spectra, we show that the optical 
QPO has a coherent twin in the X-rays. We calculate phase 
lags between the light-curves in the optical and X-rays and 
demonstrate how this information can be used to determine 
the system inclination. 

The remainder of the paper is laid out as follows. In 
Section 2 we present our data analysis and cross-correlation 
results; in Section 3 we introduce a physical model which 
can describe the simultaneous broadband variability seen in 
both X-rays and the optical and their coupled QPOs. We 
present our conclusions in Section 4. 


2 DATA ANALYSIS 
2.1 Data processing 

We reanalyse the data obtained on 2007 June 13 with 
VLT/ULTRACAM (Dhillon et al. 2007) simultaneously 
with RXTE/PCA. The data were described and published 
in Durant et al. (2009). Optical light-curves were obtained 
in two filters r' and g', with 0.143 s time resolution. Al¬ 
though the u' light-curve was also obtained, we do not use 
it because poor weather conditions during the observations 
resulted in a high noise level. X-ray data were analysed with 
the help of package heasoft. We extracted the X-ray light 
curve applying solar system barycenter correction (f axbary) 
and selected only X-rays in the energy band 3-15 keV. The 
X-ray light-curve has the same 0.143 s time resolution. 

We start our analysis of the optical light-curves with 
subtraction of the linear function fitted to the raw data. 
This procedure reduces the noise, which could affect the cal¬ 
culated characteristics, however, it does not affect principal 
results, which would have been the same even using the raw 
data. Because of a variable sky transparency during the ob¬ 
servations, we used the flux ratio of the target to the field 
star to constitute the optical light-curve. This procedure 
aims to reduce the low-frequency multiplicative atmospheric 
noise, which is the same for the target and the comparison 
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Figure 1. (a) X-ray (upper, blue), g’ (middle, green) and r' 

(lower, red) PSDs. Poisson noise has been subtracted in the X- 
ray PSD. There is a prominent optical QPO and a hint of the X- 
ray QPO is seen at the same frequency. Analysis of Durant et al. 
(2009), however, suggests that the X-ray QPO is not statistically 
significant, (b) r'/X-ray (lower) and g'/X-ray (upper, shifted up 
by 0.15 for clarity) CCFs. Orange shaded area reflects Icr errors. 
An oscillating structure is apparent at | At| > 10 s. Positive time- 
lags correspond to optical delay. 

star. It, however, does not eliminate the contamination of 
the atmospheric noise at higher frequencies corresponding 
to separation between the target and the comparison star. 
The atmospheric noise is clearly seen in the optical power 
spectral density (PSD). At some frequencies it may domi¬ 
nate over the source variability making it difficult to separate 
the intrinsic source variability from the noise contamination. 
We use arbitrary normalisation for the PSDs for the same 
reason and further consider only the PSD shape. 

2.2 Results 

We split the light-curve into 10 segments of about 290 s 
and calculate PSD in every segment. The errors are esti¬ 
mated using bootstrap method. We run 10® simulations to 
compute average PSDs from randomly-chosen 10 segments 
(each segment can be repeated), and calculate the disper¬ 
sion of these average values. In Fig. 1(a) we show optical 
and X-ray PSDs. We focus on the variability in the frequency 
range ~0.01-0.4 Hz. We note the excess in the optical power 
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Figure 2. The observed CCF (a) r'/X-ray and (b) g'/X-ray and the corresponding best-fitting cosine function. 


at frequency ~0.08 Hz (as reported in Durant et al. 2009), 
and some excess in the X-ray PSD at a similar frequency. 
Although Durant et al. (2009) state that the QPO is not 
statistically required in the X-rays. 

We calculate the CCFs in each segment by calculat¬ 
ing the inverse Fourier transform of the cross-spectral den¬ 
sity (CSD). The errors are again calculated using boot¬ 
strap method. We plot the optical/X-ray CCFs in Fig. 1(b). 
Positive time-lags correspond to optical delay. Similar to 
Durant et al. (2008), we obtain an anti-correlation of the 
two light-curves at small ( — 10 < At < 0 s) negative lags, 
accompanied by positive correlation at similarly small posi¬ 
tive lags. We note that this dip-plus-peak structure is seen 
in all the time segments, as well as in the dynamical CCF 
(figure 4 of Durant et al. 2011). In addition, we notice the 
oscillating structure at |At| > 10 s suggesting a distinct fre¬ 
quency is present in the optical and X-rays throughout the 
light-curve. Both signals are correlated at this frequency and 
the correlation holds for a number of periods. 

We fitted the CCF with the cosine function in the in¬ 
terval — 100 < At < —10 s, with frequency, phase, constant 
level and amplitude being free parameters. The results are 
plotted in Fig. 2. The best-fitting phase for both optical fil¬ 
ters was found to be zero within the chosen accuracy. The 
best-fitting oscillation frequency is f — 0.078 Hz for r' 
ray CCF, with the reduced x^/<i.o.f.=362.6/627, very close 
to the QPO frequency found by Durant et al. (2009) in the 
optical PSD. The g'/X-ray CCF shows oscillation at sim¬ 
ilar frequency of / = 0.079 Hz with xVd-o.f.=404.8/627. 
The reduced being less than unity in both cases means 
that the errors are somewhat overestimated. The fact that 
the obtained frequencies are close to the optical QPO fre¬ 
quencies gives us grounds to suspect that the QPO is also 
present in the X-ray light-curve, but is veiled there under the 
broadband noise. The significance of this result is, however, 
difficult to estimate. 

If the QPO is present in both optical and X-rays, this 
should be reflected in their combined characteristics: the 
CCF and the phase lags. We checked whether the phase-lag 
spectra demonstrate any statistically significant feature at 
the QPO frequency. This time, the total light-curve was split 
into 20 segments, about 147 s each, as we are more interested 
in having more segments to improve statistical properties 
rather than in good frequency resolution. Smooth otherwise. 



Figure 3. Phase lags between the optical and X-ray light curves 
as a function of Fourier frequency. The error bars were obtained 
using equation (9.52) of Bendat & Piersol (1986). Green squares 
and red circles with error bars show phase lags from the g' and r' 
filters compared to the X-ray light-curve, respectively. The phase- 
lags are almost independent of frequency with Acf) 2 rad at 
/ < 0.05 Hz and decrease down to A0 ~ 0.5 rad at the QPO 
frequency. 

phase lags in both filters show an abrupt drop at the QPO 
frequency by about 1.5 rad (down to A(j) ~ 0.5 rad), with 
the width of the dip of only ~0.02 Hz. To assess the signif¬ 
icance of the dip, we calculate the errors in the phase-lag 
spectra using a bootstrap method and from equation (9.52) 
of Bendat & Piersol (1986). The errors calculated by both 
methods are in good agreement in general, however some¬ 
times the former method gives larger errors due to the 27r 
phase lag uncertainty. In Fig. 3 we plot the average phase 
lags and their standard deviations calculated using the lat¬ 
ter method. We calculate the average phase lag over the fre¬ 
quency range 0.008-0.2 Hz, excluding the frequencies 0.06- 
0.1 Hz, where the QPO is likely contributing, as 

^ A(j)i/ai 

A0ave = — -, (1) 

E 1/c^i 

where A(j>i and cri are the phase-lag and its standard de¬ 
viation at i-th frequency. We obtain A())ave = 1.9 rad for 
r'/X-ray and Acfy^ve = 2.2 rad for g'/X-ray phase-lags. The 
significance of the dip is 4.7 (t for r' band and it is 3.2 ct for 
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g' band. We conclude that the dip at the QPO freqnency 
is statistically significant in both bands. Further analysis is 
applied only to r'-filter as an illustration. 

Our findings can be summarised as follows. The optical 
data of SWIFT J1753.5-0127 reveal a number of puzzles: 
there is a prominent QPO in the optical and only a hint 
of it in the X-ray PSD, a dip-plus-peak behaviour of CCF 
at small lags and its oscillating structure with the period 
close to the QPO period at larger lags, and the phase lags 
A(/) = 2 rad, almost independent of frequency with an abrupt 
drop down to nearly zero at the QPO frequency. 


3 MODELLING 

Below we develop a quantitative model to explain the fea¬ 
tures revealed by the simultaneous data. Any successful 
model should simultaneously describe the broadband noise 
both in the X-ray and the optical band, the non-trivial rela¬ 
tion between them reflected in a complex shape of the CCF 
and the phase-lags, as well as the presence of the QPOs in 
both bands with the nearly zero phase lag. The model should 
also be consistent with the known spectral characteristics of 
accreting BHs. 

There is a general agreement that the X-rays in ac¬ 
creting BHs depending on the spectral state are pro¬ 
duced by the inner hot accretion flow surrounded by a 
cold disc and/or hot non-thermal corona above the disc 
(see Poutanen 1998; Zdziarski & Gierlinski 2004; Done et al. 
2007; Gilfanov 2010; Poutanen & Veledina 2014, for re¬ 
views). There could be at least three sources of optical ra¬ 
diation: the outer irradiated accretion disc, the inner hot 
flow, and the jet (Poutanen & Veledina 2014). It is likely 
that in SWIFT J1753.5-0127 contribution of the jet is neg¬ 
ligible given a very low level of radio emission (Soleri et al. 
2010; Tomsick et al. 2015). Earlier we have developed a syn¬ 
chrotron self-Gompton-disc reprocessing model that suc¬ 
cessfully described the X-ray/optical broadband timing 
properties. In this model, the X-rays were produced by the 
inner hot flow, while the optical emission had contribution 
from two components: the synchrotron radiation from non- 
thermal electrons in the hot flow and the reprocessed emis¬ 
sion from the outer cold disc. The synchrotron component 
was assumed to be anti-correlated with the X-rays with no 
delay; such a behaviour is expected in a hot flow, when 
the increasing mass accretion rate leads to increasing syn¬ 
chrotron self-absorption in the sources. The reprocessed disc 
emission is expected to be delayed, somewhat smeared but 
positively correlated with the X-rays. 

The model described above does not consider the QPOs 
at all. Recently a physically realistic model of Lense-Thirring 
precession of the whole hot flow that describes well the X- 
ray QPOs was developed (Fragile et al. 2007; Ingram et al. 
2009). Veledina et al. (2013b) have pointed out that if the 
hot flow emits optical photons by synchrotron radiation, 
then the hot-flow precession will naturally produce not only 
X-ray but also optical QPOs, either coming in phase or anti¬ 
phase, depending on the observer position. 

The model we put forward is thus based on the com¬ 
bination of the synchrotron self-Compton-disc reprocessing 
model with the Lense-Thirring precession model of the hot 
flow. In this model, there are two distinct classes of vari¬ 


ability: the broadband noise coming from the accretion rate 
fluctuations and the (quasi-) periodic signal due to a rotat¬ 
ing emission pattern of the hot flow. We consider the hot 
flow as a processing torus, radiating both the X-rays and 
the optical emission. We see more emission at both wave¬ 
lengths when the observed surface area is maximal, resulting 
in the QPOs coming in phase. Simultaneously, the torus is 
getting brighter and fainter in an aperiodic manner depend¬ 
ing on the accretion rate and independently of the preces¬ 
sion phase. The brightening in the X-rays is accompanied 
by fading of the optical synchrotron emission as proposed 
by Veledina et al. (2013b), resulting in the anti-correlation 
of the broadband noise. The mathematical description of 
this complex variability patterns is given in the following 
sections. 


3.1 Broadband variability 

For the description of the correlated optical/X-ray broad¬ 
band variability we adopt the two-component model de¬ 
scribed in Veledina et al. (2011). The optical emission is as¬ 
sumed to be coming from two terms: from the synchrotron 
emission of the hot accretion flow which is anti-correlated 
with the X-rays and from the reprocessed emission which 
is correlated, but delayed and smeared with respect to the 
X-rays. We also now take into account the fact that the syn¬ 
chrotron emission is coming from the larger radii compared 
to the X-rays (Veledina et al. 2013a), thus it lacks the high- 
frequency signal generated at smaller radii. 

We assume that the X-ray variability is caused by prop¬ 
agating mass accretion rate fluctuations m(t) (Lyubarskii 
1997). We assume the efficient accretion and take the X- 
ray light-curve x{t) oc rh{t), thus the broadband X-ray vari¬ 
ability has the same PSD as the mass accretion rate PSD. 
We use the zero-centred Lorentzian formalism (Nowak et al. 
1999) to describe its shape. We denote the Fourier image 
of mass accretion rate by M{f), where / is the Fourier fre¬ 
quency, and use small letters with argument t to denote vari¬ 
ables in time domain and capital letters with argument / for 
those in the frequency domain hereafter. Mass accretion rate 
time-series are simulated from the observed X-ray PSD using 
the Timmer & Koenig (1995) algorithm, with zero mean. 

The optical synchrotron emission s{t) is anti-correlated 
with the X-rays, but in addition it does not have the high- 
frequency signal because it is likely produced in the outer 
part of the hot flow. This is simulated by multiplying M{f) 
by a filter function 


(///s)" + l’ 


( 2 ) 


which describes the damping of high-frequency fluctuations 
of the hot flow above the characteristic frequency fa. We thus 
define the filtered mass accretion rate mflit(t) as the inverse 
Fourier transform of the product M{f)H{f). In addition, 
the synchrotron emission can be coming earlier than the X- 
rays. This can be controlled by introducing time-lag of the 
X-ray radiation compared to synchrotron, we find, however, 
that no significant delays are required by the data, thus we 
further assume them to be zero. 

The disc reprocessed emission is simulated as a con¬ 
volution of the accretion rate and the response function 
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d{t) = m(t) r{t), where r{t) is given by 


r{t) = 


exp[-{t -tl)/t 2 ] /t 2 , 

0, t <ti. 


(3) 


It has two parameters: delay and the decay time t 2 - The 
corresponding Fonrier image of the disc emission is given by 

D{f) = Rif)M{f). (4) 


3.2 Coupling with QPO 

The QPO can be obtained from the periodic signal by 
varying either amplitnde or the oscillation frequency. Re¬ 
cent studies (Lachowicz & Done 2010) suggest that the 
QPO is composed of multiple independent oscillations of 
various duration, but the oscillation frequency is con¬ 
stant over this duration. A number of methods to cou¬ 
ple the QPO with the broadband noise were discussed 
(Burderi et al. 1997; Lazzati & Stella 1997; Menna et al. 
2003; Ingram & van der Klis 2013). Similar to the broad¬ 
band noise, we choose to simulate the QPO from its PSD 
(narrow Lorentzian) using the Timmer & Koenig (1995) al¬ 
gorithm, with zero mean. 

If the QPO is an additive process, the light-curve is a 
sum of mass accretion rate time-series and the QPO. One 
can imagine this case to be realised if a part of radiating 
matter varies only due to changing mass accretion rate and 
another part (e.g., a blob of matter) varies only due to spin¬ 
ning around a black hole. However, if the QPOs are pro¬ 
duced by the rotation of emission pattern of the entire hot 
flow (Fragile et al. 2007; Ingram et al. 2009; Veledina et al. 
2013b), the process is multiplicative, as it modulates the ac¬ 
cretion rate fluctuations, resulting in the X-ray light-curve 

a:(t) = [1-I-m(t)] [1 + £x(?(t)] - 1, (5) 

where q{t) accounts for the QPO, Ex is a positive free pa¬ 
rameter. Optical radiation is a sum of two terms 

o{t) = s{t) + rdsd{t), (6) 

where ras expresses the relative contribution of the disc and 
the synchrotron components. Power spectra of both s{t) and 
d{t) light-curves are normalised so that their integrals are 
equal, thus ras gives the ratio of the power spectra of the 
disc to that of the synchrotron. The synchrotron term is 
defined as 


s{t) = [1 - mmt(t)] [1 + £oq{t)] - 1, (7) 

where rhfiit(t) is the accretion rate light-curve, which has 
high frequencies filtered out, Eo is a positive free parameter. 
The minus sign in front of mfiit(t) accounts for the anti¬ 
correlation of broadband noise (see details in Veledina et al. 
2011), while plus signs in front of Sx and £□ aim to have 
two QPOs in phase. Thus both the X-ray and the optical 
light curves are the superposition of the broadband noise 
and the QPO variability patterns, mathematically described 
by equations (5) and (7). 

The reprocessed emission is likely modulated at the 
QPO frequency as well, however the phase shift and the os¬ 
cillation profiles require detailed knowledge of disc param¬ 
eters and the observer position (see Veledina & Poutanen 
2015). For the sake of simplicity, we consider only broadband 
fluctuations of the irradiated disc emission (equation 4). 


Table 1. Parameters of numerical modelling. 


Parameter 

Value 

^ds 

0.85 

tl 

0.1 s 

t2 

5.0 s 

fs 

0.05 Hz 

Sx 

0.9 

So 

0.8 

fo 

0.078 Hz 

A/o 

0.003 Hz 

X-ray broadband 

noise parameters^ 

A/i 

0.095 Hz 

ri 

0.13 

a/2 

1.5 Hz 

r2 

0.09 


“Adopted from Durant et al. (2009). 


3.3 Results of numerical modelling 


We start our modelling with decomposing the X-ray PSD 
into Lorentzian functions of the common form 


Hf) = 


pAfi 


n[Ap + (/ - m ' 


i = 0,l,2. 


(8) 


where Vi describe the normalisations proportional to their 
contribution to the rms. The broadband noise is described 
by two zero-centred Lorentzians, /i = /2 = 0 (after 
Nowak et al. 1999), with parameters found in Durant et al. 
(2009). The QPO is described by a Lorentzian with central 
frequency fo = /qpo found from the CCF fitting and with 
the width A/qpo found from fitting the phase-lag spectra. 
The normalisation of this Lorentzian was chosen such that 
its peak power is equal to the power of broadband variabil¬ 
ity at the same frequency, Q(/qpo) = MUqpo) (or, equiv¬ 
alently, Lo{fo) = Li(/o) + L 2 (/o)). After that, relative im¬ 
portance of the QPO is controlled by coefficients £x and Eo- 

Parameters of the model are listed in Table 1. The first 
four parameters refer to the broadband noise and affect the 
phase lags, the shape of the optical PSD and the CCF shape 
at small (|At| < 10 s) time lags, while the following four 
parameters describe the appearance of the QPO: waves in 
the CCF at large time lags, narrow spikes in the PSDs, the 
depth and the width of the dip in the phase lags. 

The disc to synchrotron ratio rus regulates the relative 
amplitude of the dip to the peak in the CCF. It simulta¬ 
neously describes how much the optical PSD is suppressed, 
relative to the X-ray PSD, at low frequencies. The latter is 
because the disc and the synchrotron emission come (almost) 
in anti-phase. Maximal suppression occurs when the syn¬ 
chrotron and disc PSDs are equal (see Veledina et al. 2011). 
The parameter also affects the phase lags: the larger is r^s 
the more phase lags resemble those of the disc (described by 
an increasing function of Fourier frequency), the smaller is 
rds the closer phase lags resemble those of the synchrotron 
term (equal to tt and independent of frequency). 

Parameters ti and t 2 define characteristic frequency of 
suppression of high frequencies in disc PSD and the fre¬ 
quency above which the disc phase lags start to substantially 
differ from zero. In the CCF, ti is responsible for the shift 
of the positive peak and t 2 determines its width. 
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Figure 4. The observed characteristics of light-curves in r' filter and X-ray (total) band. The X-ray model characteristics were obtained 
from the X-ray PSD broadband noise parameters found in Durant et al. (2009, table 5), coupled with the QPO through relation (5). 
Theoretical curves were obtained using equations (4), (6) and (7). (a) The observed CCF (red line) and the model (black line). The 
model describes well the behaviour of the CCF close to the zero lags and at negative lags. The divergence at large positive lags is likely 
caused by additional contribution of the optical QPOs coming from irradiation, which we do not account for. (b) The observed X-ray 
(blue crosses) and the optical (red crosses) PSDs and those given by the model (black solid lines). QPO is prominent in the optical and 
tentatively present in the X-ray PSD. The model well reproduces the shape of optical PSD at frequencies < 0.2 Hz. (c) The observed 
(red crosses) and the model (black line) phase lags, (d) The observed coherence 7 (/) (red bullets with error bars) as well as the arbitrary 
normalised observed (blue crosses) and the model (black line) absolute value of the CSDs. 


Parameter /s affects the width of the dip in the CCF: 
the smaller is /s the wider is the dip. We note that the dip 
shape generally depends on the shape of the filter function, 
the latter is related to the physical mechanism damping the 
high-frequency fluctuations in the hot flow. Our choice of 
H{f), however, was motivated by the fit to the CCF, rather 
than to its physical nature. 

Parameters £x and Eo simultaneously affect the QPO 
prominence in the PSDs, the wave amplitude in the CCF 
and the depth of phase lag drop at /qpo- The QPO 
Lorentzian width A/qpo not only determines the charac¬ 
teristic timescale on which the oscillations are coherent (i.e. 
at which time lags oscillations in the CCF are still seen), 
but also influences the width of the dip in the phase-lag 
spectrum which was used to determine that parameter. 

For each set of model parameters we simulate 1500 dif¬ 
ferent realisations of the optical and X-ray light curves and 
compute the average model characteristics, which we then 
fit-by-eye to the observed X-ray and optical PSDs, the CCF, 
the phase lag and the CSD (see Fig. 4). The model param¬ 
eters that describe the data well are given in Table 1. 

The model CCF consists of the correlated aperiodic 
variability and the correlated QPO, seen as oscillations at 
large negative and positive time lags (see Fig. 4a). The anti¬ 
correlation at —10 < At < 0 s is described by the syn¬ 
chrotron term, while positive peak of smaller amplitude at 


0 ;$ 20 s appears because of the disc term (see equa¬ 

tions 3, 6 and 7). Interplay of the two optical terms leads to 
a sharp feature at At = 0, where their relative role is rapidly 
changing. The CCF shape is not well described at larger pos¬ 
itive lags, possibly because of the additional presence of the 
optical QPO from the disc, which we do not account for. A 
significant excess at the QPO frequency is seen is also seen 
in the CSD (Fig. 4d), the Fourier image of the CCF, sup¬ 
porting our finding that the oscillations at this frequency in 
the optical and the X-rays are correlated. 

The X-ray PSD is obviously well described by the model 
just because we used the parameters from Durant et al. 
(2009) that fit it (see Fig. 4b). The optical PSD shape is 
consistent with the model at frequencies / < 0.2 Hz. We 
note that while the X-ray PSD is almost constant at fre¬ 
quencies / < 0.03 Hz, the optical power is somewhat sup¬ 
pressed, both in the model and in the data. The suppression 
again comes from the interplay of two optical components 
which nearly cancel each other at sufficiently low frequencies 
as they come in anti-phase. 

The model phase lags are close to 2 rad at all frequencies 
except of the QPO (see Fig. 4c). They reflect a joint con¬ 
tribution of the synchrotron {A(j) = n rad) and of the disc 
components (Ac^ ~ 0 in the frequency range of interest) to 
the phase lags. The dip in phase lags at the QPO frequency 
is also well reproduced. The phase lags do not reach zero in 
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Figure 5. (a) The PSDs of SWIFT J1753.5—0127 (black squares), of the comparison star (cyan circles) and of the ratio of their fluxes 
(red crosses, overplotted on both PSDs) for the night 2007 June 13. The shape of the ratio-PSD matches the PSD of the comparison 
star at frequencies / > 0.2 Hz, but is substantially different from the PSD of the source alone, (b) The PSD of SWIFT J1753.5-0127 for 
night 2007 June 18 (blue triangles), as compared to the PSD of the target/comparison star ratio on night 2007 June 13 (red crosses). 
The normalisation is chosen such that the two PSDs match at the QPO frequency. They also match up to frequency ~ 0.2 Hz and 
substantially deviate at higher frequencies. 


spite of the fact that the QPOs are in phase, because there is 
still a substantial contribution of the broadband variability 
with A0 ~ 2 rad even at the QPO frequency. 

The QPO is apparent in the model of optical PSD and 
there is a hint of the QPO in the model of X-ray PSD, con¬ 
sistent with the data. The model significantly underpredicts 
the optical PSD above frequencies ~ 0.2 Hz. The cut-off in 
the model PSD is related to the width of the anti-correlation 
in the CCF, requiring the synchrotron term to have high fre¬ 
quencies damped. The detailed shape and width of the dip 
in the CCF depends on the choice of the filter function: the 
smaller is the characteristic filter frequency the wider is the 
dip and the less power is stored at high frequencies. At the 
same time, the disc term does not show much power at high 
frequencies either. 

Some of the model parameters can be determined from 
the data with high accuracy. For example, the QPO fre¬ 
quency and its width are well defined. On the other hand, 
is rather uncertain. We fit the data with t\ = 0.1 s, while 
an equally good fit can be achieved with up to 0.5 s. 
The difference can be only seen at high frequencies, which 
are likely contaminated by the noise thus are not reliable to 
prove parameter. Although most of the parameters do not 
have direct physical interpretation, parameters t\ and t 2 are 
likely related to the characteristic size of the reprocessing 
region. Interestingly, the value of f 2 = 5 s is consistent with 
the maximum time delays corresponding to the light travel 
back and forth the accretion disc of radius 6 x 10^° cm ex¬ 
pected in this system (Neustroev et al. 2014). 


3.4 Nature of the optical variability above 0.1 Hz 

Although the presented model describes rather well the X- 
ray and the optical PSD, the CCF, the phase lags and the 
CSD, there are some deviation in the optical PSD at high fre¬ 
quencies. The nature of the high-frequency variability above 
> 0.1 Hz is not certain. Because the model reproduces the 
shape of the entire CCF (which reflects the correlated sig¬ 
nal) and the CSD, the variability is unlikely intrinsic to the 


source. The high-frequency noise might be due to the atmo¬ 
sphere, not completely removed by taking the ratio of the 
BH to the comparison star fluxes. 

We compare the target and comparison star r'-filter 
PSDs with the PSD of target/comparison star ratio in 
Fig. 5(a). We note that the PSD of the ratio is almost iden¬ 
tical to the PSD of the comparison star (which reflects the 
atmospheric noise fluctuations) at / > 0.3 Hz. At the same 
time, the shape of the target PSD is substantially different, 
demonstrating a hump at these frequencies. In Fig. 5(b) we 
compare PSDs of the target and the ratio taken on 2007 
June 13 with the target PSD taken on 2007 June 18, when 
the weather conditions were much better. The PSDs are nor¬ 
malised so that they match at the QPO frequency. We see 
that the June 18 target and June 13 ratio PSDs match ev¬ 
erywhere below / < 0.2 Hz. Assuming the target PSD did 
not suffer substantial change within these days, we expect 
the shape of target June 13 PSD to be different from the 
one demonstrated by the target/comparison star ratio PSD. 
We also note drop of the coherence towards high frequencies 
(see Fig. 4d), meaning the decrease of power of the correlated 
signal in the PSDs. These facts support our suggestion that 
the high-frequency optical power is due to the atmospheric 
noise rather than being intrinsic to the source. 

We note that this noise may also “veil” the high- 
frequency part of the source optical PSD. This affects the 
CCF shape making both the dip and the peak broader, 
as the high-frequency correlated signal is lost, and can be 
the reason for a much narrower dip and peak detected in 
the 2008 simultaneous data (Durant et al. 2011). From the 
model point of view, the synchrotron and the disc compo¬ 
nents may vary on shorter time-scales than those obtained 
from modelling, thus no filter function is needed. Detailed 
comparison of the data taken in 2007 and 2008 is needed to 
clarify this question. 

Alternatively, the high-frequency component might be 
related to synchrotron photons from the hot flow, which are 
not efficiently Comptonized (e.g., from the outer regions of 
the hot flow which are more transparent to the synchrotron), 
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and thus are not correlated with the X-rays. Given a very 
low level of radio emission during the outburst history of 
the source (Soleri et al. 2010; Tomsick et al. 2015), any jet 
contribution to the optical wavelengths is doubtful. 


4 CONCLUSIONS AND OUTLOOK 

We presented the analysis of the simultaneous optical/X-ray 
data from BH binary SWIFT J1753.5-0127. We found pres¬ 
ence of the QPO in the X-ray data and demonstrated that it 
is coupled with the reported earlier optical QPO at the same 
frequency. We showed that the oscillations are coherent, as 
the periodic signal is seen in the CCF. The best-htting os¬ 
cillation frequency obtained by approximating the CCF at 
large lags with the cosine function gives /qpo = 0.078 Hz, 
the same frequency as the QPO found in the optical PSD. 
We found the phase lags are almost independent of the 
Fourier frequency, A(j) ~ 2 rad, but demonstrated a signif¬ 
icant (4.7(7 and 3.2a for r'/X-ray and p'/X-ray phase lags, 
respectively) decrease at the QPO frequency. This suggests a 
separate narrow component is present both in optical and X- 
ray PSDs, with the phase lags significantly lower than those 
of the broadband noise. The fact that the CCF was well fit¬ 
ted with the cosine function with zero phase suggests that 
these narrow components are coming in phase {A<f) = 0). 
These hndings not only extend the range of X-ray luminosi¬ 
ties where the QPOs are present, suggesting those might 
be hidden in the noise during the hard state, but also put 
strong constraints on the origin of optical QPOs, which has 
been disputed since their very discovery. 

We developed a model which simultaneously explains 
the shape of the CCF, phase lags and PSDs. It is based on 
the temporal behaviour of the hot accretion flow proposed 
in Veledina et al. (2011, 2013b). In this model, the aperiodic 
variability in the X-rays is assumed to reflect the propagat¬ 
ing fluctuations in mass accretion rate. The optical emission 
is composed of two terms: the synchrotron component which 
varies in anti-phase with the X-rays and lacks high-frequency 
signal (as it is coming from larger radii), and the irradiated 
disc component which is delayed and smeared with respect 
to the X-rays. Their interplay results in phase-lags of 2rad, 
almost independent of frequency. The X-ray and optical syn¬ 
chrotron broadband variability are coupled to the QPO in 
a multiplicative way, as in the model of the Lense-Thirring 
precession of the hot flow. The low value of the phase lags 
at the QPO frequency is caused by the superposition of the 
~2 rad lag coming from the aperiodic variability and the 
QPOs which are coming in phase, A((> = 0. 

The fact that the QPOs in the optical and X-rays come 
in phase restricts in the model of Veledina et al. (2013b, see 
their fig. 5) the orbital parameters of the system, suggesting 
inclination i > 50°. This is consistent with the constraint 
i > 40° obtained from the prominent variability of the optical 
flux and of the line equivalent width at the orbital period 
(Neustroev et al. 2014). If indeed i > 50°, for the measured 
mass function f{M) < 0.95Mq (Neustroev et al. 2014), this 
implies the BH mass of < 2.5M0, strongly restricting its 
possible progenitors. 
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